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ABSTRACT 

This paper presents performance results for pulse detonation engines taking into account the effects of dissocia- 
tion and recombination. The amount of sensible heat recovered through recombination in the PDE chamber and 
exhaust process was found to be significant. These results have an impact on the specific thrust, impulse and fuel 
consumption of the PDE. 


INTRODUCTION 

In a previous publication (Ref. 1) it was shown that the high temperatures associated with detonation caused a 
substantial amount of dissociation, yielding high concentrations of intermediate species. This process was 
accompanied by an 11% decrease in the amount of energy relative to a ramjet available for the production of thrust. 
The consequence of this loss was to reduce the thermal efficiency of the PDE below the value of the Brayton cycle 
efficiency at low values of ram air temperature ratios (Ref. 1). Although recombination effects were qualitatively 
considered, they were not evaluated previously. In this paper, a quantitative evaluation of the effect of 
recombination was carried out and the amount of sensible heat recovery was determined. In this manner, it may be 
determined if any significant changes are required to the previous calculations of energy release. The work in 
references 1 and 2 was based on chemical equilibrium whereas in the present study, finite rate calculations were 
employed. 


PHYSICAL MODEL AND NUMERICAL METHOD 


Geometry and Test Gas 

Computations were performed for a constant area detonation tube having a length of 1.0 meter and a diameter of 
6.6 centimeters. The tube is closed on the left and open on the right. Hydrogen and air having an equivalence ratio, 
4>, of 1.0 and 0.6, with a temperature of 298 K were used as the reactants. Calculations were performed for zero 
forward flight velocity and at a pressure of 1 bar. At the start of the computation, the entire tube is filled with the 
hydrogen-air mixture. The ambient pressure outside the detonation tube is set at 1 bar. 

To initiate the detonation, we use a small region (2 cm long from the detonation tube head end) of high tempera- 
ture and high pressure (2500 K and 40 bar) nitrogen gas. In this direct initiation method, the shock wave generated 
at the interface will transition into a Chapman-Jouguet detonation wave as shown later. 


* Chief Scientist, Turbomachinery and Propulsion Systems Division, 
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Computational Method 

Finite rate calculations were used to compute the species evolution and resulting thrust and impulse. The 
analysis was carried out using an in-house developed time-accurate CFD code (Refs. 4, 5). The code solves the 
axisymmetric Navier-Stokes equations including finite-rate chemistry and real gas effects using an implicit, total 
variation diminishing (TVD) algorithm. It includes a generalized detailed chemistry capability, various options for 
turbulence models, and steady-state or time accurate marching algorithms. In the present study, viscous effects were 
not included, and therefore, the Euler equations were actually used. 

The numerical method used for solving the governing equations is described in detail in Ref. 4, and briefly sum- 
marized here. The equation set is discretized using a fully implicit, first-order-accurate in time, variable-step back- 
ward differentiation formula (BDF) method. The numerical fluxes are discretized using a second-order spatially 
accurate TVD scheme. The resulting equations are then linearized in a conservative manner and solved iteratively, 
by using a lower-upper relaxation procedure consisting of successive Gauss-Seidel (LU-SGS) sweeps. 

The chemical reaction mechanism for hydrogen-air combustion was based on Jachimowski’s model (Ref. 6), 
and consists of 27 elementary reactions among 12 species. The inversion of large matrices is avoided by partitioning 
the system into reacting and nonreacting parts. Consequently, the matrices that have to be inverted are of the same 
size (N-N, where N is the number of reacting species) as those that arise in the commonly used point implicit 
methods. An important advantage of the present method is that, because it is fully implicit, it is stable for large 
values of the CFL number, thereby enabling the use of relatively large time steps to minimize computational cost. 

In order to maintain good resolution of the detonation front at all times, without having to use thousands of grid 
points, a multi-level, dynamically adaptive grid was implemented in which a very fine subgrid continuously slides 
along with the detonation wave front (Ref. 7). The minimum and maximum grid spacings used in the present study 
were 0.00195 cm and 1.0 cm, respectively. 


RESULTS 

Finite rate chemistry computations were performed for (f)=1.0 and 0.6. The resulting pressure contours are shown 
in figure 1 for the stoichiometric case. Figures 1 through 8 show the (f)= 1 .0 results. The initial conditions are shown 
in Fig. 1(a), where the small high pressure zone is observed at the closed end of the tube. At the interface of this 
region and the combustible mixture, an expansion wave propagates toward the closed end and is reflected from it. 
Simultaneously, a shock wave moves towards the open end of the tube compressing and heating the hydrogen-air 
mixture. After a short induction time, combustion begins at the interface boundary sending compression waves 
upstream and downstream. The compression waves overtake the shock, accelerating it. The shock and combustion 
front become coupled, forming a detonation wave. Figure 1(b) shows the detonation wave near the end of the tube. 
The next plots (Figs. l(c)-(f)) show the detonation expansion process into the surrounding air. 

It may be seen that positive pressure is present for approximately 3 milliseconds, during which thrust can occur. 
The arrival of the rarefaction wave at the head end of the tube can be seen occurring between the last two time 
sequences. 

The computed detonation velocity as a function of time is shown in figure 2, and compared with predictions 
from the CEA equilibrium code of Gordon and McBride (Ref. 8) for a Chapman-Jouguet detonation. It is observed 
in figure 2 that after an initial overshoot during the short transient phase, the detonation speed reaches a nearly 
constant value that is 0.8% higher than that predicted by the CEA code. The corresponding sensible heat release is 
shown in figure 3 for both the finite rate calculations and the chemical equilibrium results from the CEA code. For 
comparison, the heat of reaction for hydrogen is also plotted. Note that the CFD results represent the heat released 
per mass of fuel burnt. It may be seen in figure 3 that the finite rate and equilibrium results differ by 7.5%. The 
higher heat release for the finite rate calculation is due to chemical recombination inside the tube. Note that both the 
detonation speed and the heat release are increasing slowly for t > 0.05 msec. Therefore, the higher detonation speed 
computed with the CFD code, as compared with the equilibrium value, could be a result of finite rate chemistry 
effects. 
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Figure 4 shows that the passage of the initial detonation pressure spike and temperature rise occurred at 
0.5 msec, followed by a rapid decay within 1 msec. The plateau region, where the pressure remains at a level value 
is seen to persist up to about 3 msec as indicated previously in figure 1 . These properties are shown at the exit of the 
detonation tube. 

The changes in gas species are shown in figure 5. The species concentrations are seen to decay rapidly within 

2 msec, followed by a slowly decreasing trend. The concentrations of OH, H02, and H202 decline by an order of 
magnitude, whereas the O and H decline about 2 orders of magnitude in the first 2 msec. 

In a second set of computations, the finite rate results were compared with the case where the flow in the tube 
was frozen when the detonation wave reached the exit of the tube. Figure 6 shows that the temperature profile for 
the frozen case is approximately 7% lower than that obtained with the finite rate. The species mole fractions at the 
exit of the tube are shown in figures 7(a)-(c). It is seen that the species concentrations of the finite rate results are 
lower than those occurring with the frozen flow assumption with the exception of NO. 

The corresponding thrust and impulse for the computations described above are shown in figure 8. Initially, the 
sharp rise in pressure creates a short duration thrust spike (Fig. 8(a)), which is quickly followed by a longer plateau 
region of 2 msec duration and a subsequent decay to zero. It is primarily during this level pressure duration that 
PDE thrust is generated. The finite rate profile is slightly greater than that of the frozen case from the range of 1 to 

3 msec. The corresponding values of impulse are shown in figure 8(b). The value of the impulse at its maximum 
point is greater by 5% from the corresponding frozen case. 

The PDE results for an equivalence ratio of 0.6 are shown in figures 9 through 12. Figure 9 shows a comparison 
of the sensible heat release for the finite rate calculation and the chemical equilibrium results for both equivalence 
ratios. It is noted that the lean equivalence ratio results indicate that the dissociation and recombination processes 
are less important on their effect on sensible heat release. 

Figure 10 shows the temperature profile as observed at the detonation tube exit location. The temperature for 
both frozen and finite rate chemistry show only a minor difference. A comparison of the intermediate species are 
shown in figure 11. Generally, the concentrations of the species that contribute to dissociation losses are 
significantly lower for cf)=0.6 than those for the stoichiometric results. The corresponding force and impulse 
behavior are shown in figure 12. In contrast to the cf)= 1 .0 results, no substantial differences are present. 

The fuel specific impulse was also determined and is shown in figure 13(a) for both equivalence ratios. The 
specific impulse for cf)=0.6 was higher than that for cf)= 1 .0. However, the thrust for cf)=0.6 was lower than that at 
(])=1.0 (see figures 12(a) and 8(a)). In figure 13(b), the starting transient that occurs in the first 0.05 msec was 
neglected. This is justified on the basis of the need to introduce an ignition mechanism to obtain detonation which 
creates an artificial force that is not representative of the actual initiation process. The specific impulse result based 
on this correction is shown in figure 13(b). It is observed that the slopes for each equivalence remains unchanged, 
whereas the maximum values decreased when the ignition transient was neglected. The impulse level also decreased 
by approximately 10% in the absence of the transient. As in previous results, figure 13(a), the cf)=0.6 results show 
greater specific impulse than that for (])=1.0. 


CONCLUSIONS 

The results show that the assumption of frozen flow subsequent to the detonation wave yields lower 
temperatures at the tube exit for cf)= 1 .0; whereas in the case of cf)=0.6, the differences are negligible. Significant 
differences were observed in the species mole fractions depending on the equivalence ratio. In general, the rapid 
decay within the first 2 msec was followed by a more gradual decline for both equivalence ratios. As expected, the 
absolute level of dissociation was higher for the stoichiometric case. In addition, the force exerted at the head end of 
the detonation tube shows a small difference, which occurs during most of the process for the cf)= 1 .0 case. This force 
difference was much smaller for the cf)=0.6 case. The impact on the maximum impulse is of the order of 5.0% for the 
4>= 1 .0 case, and 2.8% for the (])=0.6 case. 
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The main conclusion drawn from this study was that the dissociation losses occurring during the detonation pro- 
cess for cf)= 1 .0 reduces the sensible heat release by approximately 16.7% relative to the heating value, when 
chemical equilibrium was assumed. Recombination in the PDE reduces this loss to 10.8%. For 4>=0.6, dissociation 
effects reduced the sensible heat release by 4.7% (when chemical equilibrium was assumed), and recombination 
reduces the loss to only 0.6%. 

The equivalence ratio comparison indicates that the PDE performance is improved by operation at off- 
stoichiometric conditions. This suggests operation of a PDE on the lean side provided that sufficient thrust is 
produced. 
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Figure 1. Pressure contours showing temporal evolution of the detonation wave. H ,-air mixture, (f)= 1 .0, 
p 0 =1.0 bar, T o =298 K. 
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Figure 2. Detonation speed variation with time. H,-air mixture, 4>= 1 .0, p 0 =1.0 bar, T o =298 K. 



Figure 3. Heat release per mass of burnt fuel as a function of time, H,-air mixture, cf)= 1 .0, p o =1.0 bar, T o =298 K. 
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Figure 4. Flowfield conditions at the exit of the detonation tube (x = 1.0 m) as a function of time ( <f)= 1.0). 



Figure 5. Species mole fractions at the exit of the detonation tube (x = 1.0 m) as a function of time (<|>=1.0) 
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Figure 6. Comparison of temperature at the exit of the detonation tube (x = 1.0 m) between finite rate chemistry 
flow, and flow that has been frozen after the detonation wave exits the tube (4>= 1 .0). 
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Figure 7. Comparison of species concentration at the exit of the detonation tube (x = 1.0 m) between finite rate 
chemistry flow, and flow that has been frozen after the detonation wave exits the tube (<|>=1.0). 
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Figure 8. Comparison of force (a) and impulse (b) at the exit of the detonation tube (x = 1.0 m) between finite rate 
chemistry flow, and flow that has been frozen after the detonation wave exits the tube (cf)= 1 .0). 
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Figure 9. Heat release per mass of burnt fuel as a function of time for the two equivalence ratios. 
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Figure 10. Comparison of temperature at the exit of the detonation tube (x = 1.0 m) between finite rate chemistry 
flow, and flow that has been frozen after the detonation wave exits the tube (cf)=0.6). 
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Figure 1 1 . Comparison of species concentration at the exit of the detonation tube (x = 1.0 m) between the finite rate 
chemistry computations for <))= 1 .0 and cf)=0.6. 
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Figure 12. Comparison of force (a) and impulse (b) at the exit of the detonation tube (x = 1.0 m) between finite rate 
chemistry flow, and flow that has been frozen after the detonation wave exits the tube (cf)=0.6). 
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Figure 13. Comparison of specific impulse between the finite rate chemistry computations for cf)= 1 .0 and 4>=0.6. In 
part (b), the force produced during the starting transient phase (first 50 jisec) was neglected. 
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